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Abstract 

The moss Physcomitrella patens is part of an early divergent clade of land plants utilizing the plant hormone cytokinin 
for growth control. The rate-limiting step of cytokinin biosynthesis is mediated by isopentenyltransferases (IPTs), 
found in land plants either as adenylate-IPTs or as tRNA-IPTs. Although a dominant part of cytokinins in flowering 
plants are synthesized by adenylate-IPTs, the Physcomitrella genome only encodes homologues of tRNA-IPTs. This 
study therefore looked into the question of whether cytokinins in moss derive from tRNA exclusively. Targeted gene 
knockout o1 ipti (d|ipt1) along with localization studies revealed that the chloroplast-bound IPT1 was almost exclu- 
sively responsible for the A37 prenylation of tRNA in Physcomitrella. Ultra-performance liquid chromatography-tandem 
mass spectrometry (UPLC-MS/MS)-based cytokinin profiling demonstrated that the total amount of all free cytokinins 
in tissue was almost unaffected. However, the knockout plants showed increased levels of the A/6-isopentenyladenine 
(iP)- and frans-zeatin (tZ)-type cytokinins, considered to provide active forms, while c/s-zeatin (cZ)-type cytokinins 
were reduced. The data provide evidence for an additional and unexpected tRNA-independent cytokinin biosynthetic 
pathway in moss. Comprehensive phylogenetic analysis indicates a diversification of tRNA-IPT-like genes in bryo- 
phytes probably related to additional functions. 
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Introduction 

Cytokinins (Cks) are A'^-substituted adenine derivatives act- which can be divided into an aromatic and an isoprenoid 

ing as phytohormones in various developmental processes in group depending on the side chain coupled to the of the 

plants. Natural Cks represent a large group of phytohormones adenine (Strnad, 1997; Mok and Mok, 2001). The isoprenoid 



Abbreviations: Ck, cytokinin; cZ, c/s-zeatin; cZOG, c;s-zeatin 0-glucoside; cZR, c/s-zeatin riboside; cZRS'MP, c/s-zeatin riboside-5'-monoptiosphate; cZROG, 
c/s-zeatin riboside 0-glucoside; DHZ, diliydrozeatin; DHZOG, ditiydrozeatin 0-glucoside; DHZR, dihydrozeatin riboside; DHZRS'MR ditiydrozeatin riboside- 
5'-monophosptiate; DHZROG, dihydrozeatin riboside 0-glucoside; DMAPR dimethylallyl diptiosphate; HMBDR hydroxymettiylbutenyl diphosptiate; IP, 
/V^-isopentenyladenine; iP9G, /V^-isopentenyladenine-9-glucoside; iPR, /V^-isopentenyladenosine; IPRS'MP, /V^-isopentenyladenosine-S'-monophosphate; tZ, frans- 
zeatin; tZOG, frans-zeatin 0-gluooside; tZR, frans-zeatin riboside; tZR5'MP, frans-zeatin riboside-5'-monoptiosphate; tZROG, frans-zeatin riboside O-glucoside. 
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Cks occur as four types, the A^-isopentenyladenine (iP), cis- 
zeatin (cZ), trans-zeatin (tZ), and dihydrozeatin (DHZ) types, 
varying in abundance as well as stereochemistry of the side 
chain hydroxyl group and/or in the saturation of the side 
chain in the case of DHZ types. Furthermore, Cks can be 
vastly modified, either at the A^^, A^^, or position of the 
purine ring or at the terminal side chain hydroxyl group. 
Major modifications are glycosylation, aminoacylation, and 
phosphorylation (Mok and Mok, 2001). The different types 
and conjugates of Cks can differ greatly in their biological 
activity and abundance (Sakakibara, 2006). 

Cytokinin biosynthesis 

The rate-limiting step in biosynthesis of isoprenoid Cks in 
plants is catalysed by isopentenyltransferases (IPTs), add- 
ing an isoprenoid side chain to the A^-amino group of an 
adenine nucleotide. Depending on the nature of the nucle- 
otide substrate, two different forms of IPTs are known in 
plants (Fig. 1) (Kamada-Nobusada and Sakakibara, 2009). 
Adenylate-IPTs use either ATP, ADP, or AMP as a substrate 
and are encoded in multiple copies in all flowering plants 
(Kakimoto, 2001; Takei et al., 2001). The adenylate-IPTs 
of Arahidopsis are mainly responsible for the synthesis of 
iP- and tZ-type Cks (Miyawaki et al, 2006). Interestingly, 
the first adenylate-IPT to be described was cloned from the 
plant pathogen Agrobacterium tumefaciens and was shown to 
be the factor driving Agrobacterium-inAviCQd tumorigenesis 
in plant root tissues (Barry et al, 1984; Lichtenstein et al, 
1984), and proved to be functional in the moss Physcomitrella 
patens (Reutter et al , 1998). 

tRNA-IPTs (tRNA delta2-isopentenylpyrophosphate 
transferases; also known as IPPTs; EC 2.5.1.75) on the other 
hand can be found in all organisms except Archea. They cata- 
lyse the addition of dimethylallyl diphosphate (DMAPP) 
to a tRNA-bound adenine nucleotide, 3' adjacent (A37) to 
the anticodon of UNN-decoding tRNAs (Murai, 1994; 
Miyawaki et al, 2006). This A37 modification increases the 
precision of protein biosynthesis by stabilizing codon-anti- 
codon binding (Konevega et al, 2006). In plants, the tRNA- 
bound biosynthetic pathway contributes to the formation of 
certain Cks by releasing them during tRNA breakdown. In 



Arahidopsis, tRNA-IPTs are exclusively responsible for the 
synthesis of cZ-type Cks (Miyawaki et al, 2006). 

Cytol<inins in Physcomitrella 

The last common ancestor (LCA) of mosses and flowering 
plants lived ~500 million years ago (MYA) (Lang et al. , 2010). 
Consequently, development of the moss Physcomitrella 
is mainly controlled by the evolutionarily old phytohor- 
mones auxin, Ck, and abscisic acid (for reviews, see Decker 
et al, 2006; Cho et al., 2009; von Schwartzenberg, 2009). 
Physcomitrella has been established as a model system in evo- 
devo research and plant molecular genetics due to the avail- 
ability of a fully sequenced genome as well as exceptionally 
high rates of homologous recombination enabling reverse 
genetics (Schaefer and Zryd, 1997; Kamisugi et al, 2006; 
Lang et al, 2008; Rensing et al, 2008; Prigge and Bezanilla, 
2010). Thus the evolution of phytohormone signalling and 
metabolism can be studied in comparative analyses utiliz- 
ing moss (Vandenbussche et al, 2007; Paponov et al, 2009; 
Kopecna et al, 2013). 

Physcomitrella is the first bryophyte with a reported Ck 
profile (von Schwartzenberg et al, 2007). In protonema 
from liquid cultures, the cZ-type cytokinins are dominant, 
with cZROG (c«-zeatin riboside 0-glucoside) accounting 
for 80% of the overall Ck pool. The other Ck types, iP, tZ, 
and DHZ, were also detectable at intracellular concentra- 
tions in the picomolar range. In the culture medium, the iP 
conjugates were found to be the most abundant type, with 
7V*-isopentenyladenosine-5'-monophosphate (iPR5'MP) as a 
main constituent. 

Cks are of particular importance for the progression of the 
life cycle of mosses as they induce the production of three- 
faced apical cells (buds), which later develop into the leafy 
shoots (Cove and Ashton, 1984). An external appHcation of 
Cks results in a strong overproduction of buds (Reski and 
Abel, 1985), and internal Ck accumulation precedes bud for- 
mation (Schulz et al, 2000). However, not all Cks exert the 
same hormonal activity when applied to protonema cultures. 
In budding bioassays, the free bases iP, tZ, and benzyladenine 
(BA) are the most effective, while their corresponding ribo- 
sides are less active. Strikingly, in the same assays, none of 




AMP/ADP/ATP DMAPP/HMBDP 

Fig. 1. Simplified sclieme of ihe two alternative cytokinin (Ck) biosynthetic pathways in plants, (a) Prenylatlon of UNN-decodlng tRNAs at adenine (A37) by 
tRNA-IPTs and subsequent release of cytokinin nucleotides by tRNA degradation, (b) Direct synthesis by adenylate-IPTs. DMAPP, dimethylallyl diphosphate; 
HMBDP, hydroxymethylbutenyl diphosphate. Data taken from Sakakibara (2006). (This figure Is available In colour at JXB online.) 
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the tested cZ conjugates had a bud-inducing effect, leading to 
the conclusion that cZ-type Cks are not contributing to bud 
induction (von Schwartzenberg et al, 2007). This apparent 
lack of hormonal activity stands in contrast to the dominant 
abundance of cZ-type Cks in Physcomitrella tissue. 

The full genomic sequence of Physcomitrella (Rensing 
et al, 2008; Zimmer et al., 2013) reveals seven ipt genes, 
which are at first sight all to be considered as tRNA-IPTs 
based on their sequence homology to other plant tRNA-IPTs 
(Sakakibara, 2006; Frebort et al, 2011). IPTl has been stud- 
ied by Yevdakova and von Schwartzenberg (2007). A yeast 
complementation assay confirmed the IPTl protein to be a 
functional tRNA-IPT. Findings such as the dominant occur- 
rence of cZ-type Cks in tRNA and whole tissue extracts of 
Physcomitrella, together with the sequence homology of all 
IPTs to known, functional tRNA-IPTs, led to the working 
hypothesis that the tRNA biosynthetic pathway is of great 
importance in Physcomitrella (Yevdakova et al., 2008; von 
Schwartzenberg, 2009). The remaining six IPTs have only 
rarely been studied, but Patil and Nicander (2013) have 
shown recently that IPT4 and IPTS also possess tRNA-IPT 
function. 

In this study, insights are provided into the origins of Cks 
in Physcomitrella focusing on IPTl by the generation of tar- 
geted gene knockout plants. The results of the phenotypic 
analysis, detailed analysis of free and tRNA-bound Cks, and 
a phylogeny of IPT proteins reveal unexpected, multiple ori- 
gins of Cks and indicate that the in vivo function of anno- 
tated tRNA-IPTs has to be re-assessed in Physcomitrella and 
possibly other organisms. 



Materials and methods 

Plant materia! and growth conditions 

The Physcomitrella patens (Hedw.) Bruch & Scliimp wild type 
used in this study was originally collected from Grandsden Wood, 
Huntingdonshire, UK (1968) and is the clone that was sequenced. 
Under standard growth conditions, strains were kept at 25 °C, in 
white light (Phillips TML, Hamburg, Germany) at 100 |iE m"^ s"' 
under a hghtidark cycle of 16:8 h. Tissue of Physcomitrella grown in 
liquid culture was used for transformation of protoplasts, Ck pro- 
filing, and budding assays. Liquid A'BCD(N)TV media [0.356mM 
Ca(N03)2, 1.01 mM MgS04, 1.84mM KH2PO4, lOmM KNO3, 
0.044mM FeSOJ supplemented with Hoagland trace element 
solution (1ml F') and the vitamins nicotinic acid (8 \iM), /i-amin- 
obenzoic acid (1.8 nM), and thiamine HCl (1.5 |j,M) were prepared 
according to Wang et al. (1981). For growth assays, Physcomitrella 
was grown in Petri dishes on ABCTV agar medium according to 
Knight (1988). 

Vector construction of IPT1::GFP and localization by confocal 
laser scanning microscopy (CLSM) 

The gene models of the moss IPT family were manually curated 
based on cDNA and expressed sequence tag (EST) evidence using the 
genome browser and annotation service of the moss model organism 
database (http://www.cosmoss.org; Lang et al, 2005). The complete 
open reading frame of PpIPTl (GenBank accession no. EF512463.1; 
curated cosmoss.org gene model Ppls96_115V6 lindner.l; exclud- 
ing the stop codon) was amplified from pFL61_PpIPTl (Yevdakova 
and von Schwartzenberg, 2007) with the primers iptl_for(^vrII), 



CGCGCCTAGGATGGTGAGTTTGCAGTTTGAG; and iptl_ 
rev(5flmHI), ATGCGGATCCACATAGAGGTCGTCAAGGATG, 
digested with Avrll and BamHl, and cloned into the pLNU vec- 
tor (http://www.dna-cloning.com/vectors/Vectors_with_markers/ 
pLNU-GFP.gb) 5' of the green fluorescent protein (GFP)-encoding 
sequence (the pLNU-GFP_Ppiptl vector card is given in 
Supplementary Fig. SI available at JXB online). The plasmid was 
then transfected into Physcomitrella protoplasts by polyethylene 
glycol (PEG)-mediated transfer according to Schaefer et al (1991). 
GFP fluorescence in transiently transformed protoplasts was moni- 
tored 5 d after transfection using a TCS SPE confocal laser scanning 
microscope, and images were acquired with the LAS AF lite soft- 
ware (both Leica, Wetzlar, Germany). 



Generation and characterization ofd\lpt1 mutants 

IPTl -deficient mutants were generated by gene targeting. The 
gene-disrupting vector pBNR_PpIPTl, based on pBNRr (Schaefer 
et al, 2010), carries two ~1000bp 5' and 3' genomic fragments of 
the native iptl locus (deletion of the Ppl assembly region scaf- 
fold_96:726668.. 728288), flanking a 35S:nptII resistance-mediating 
cassette. Recombinant transformants were obtained in three cycles 
of selection on ABCNTV agar medium containing G418 and were 
characterized by PCR and reverse transcription-PCR (RT-PCR) 
for iptl gene disruption as well as loss of transcript (Supplementary 
Fig. S3 at JXB online). The selected mutant moss lines are named 
d|iptl#9, d|iptl#10, and d|iptl#13, and have been deposited at the 
International Moss Stock Center (http://www.moss-stock-center. 
org/) with the respective IMSC accession numbers 40697, 40698, 
and 40699. 



Cytokinin measurements in tissue, medium, and tRNA samples 
Physcomitrella wild type and three independent d|iptl mutants were 
grown in liquid cultures under standard conditions as previously 
described in von Schwartzenberg el al (2007). Prior to prepara- 
tion of cultures, the tissue used for inoculation was washed with an 
excess of fresh medium. For sampling, 200 ml aliquots were taken 
after 22 d of culture. Tissue and medium from the same sample were 
separated by filtration. The tissue was immediately frozen in liquid 
nitrogen, and medium and tissue were stored at -20 °C for a maxi- 
mum of 2 weeks prior to freeze-drying. 

For analysis of endogenous Cks, extraction and purification were 
performed according to the method described by Novak et al. (2003). 
The samples were purified using a combined cation (SCX-cartridge), 
anion (DEAE-Sephadex-C18-cartridge) exchanger and immunoaf- 
finity chromatography (lAC) based on a wide range of specific mon- 
oclonal antibodies against Cks. The levels of Cks were quantified 
by ultra-performance liquid chromatography-electrospray tandem 
mass spectrometry (UPLC-MS/MS) using stable isotope-labelled 
internal standards as a reference (Novak et al, 2008). 

In order to obtain sufficient quantities of tissue for Ck determi- 
nation in tRNA, cultures independent from those used for free Ck 
profiling were grown under the same conditions for 22 d. 

Extraction and purification of tRNA was performed according 
to a protocol described by Maass and Klambt (1981), including 
modifications described by Stirk et al (201 1). Aliquots of 60 |.i,g of 
tRNA were hydrolysed with 2M KOH overnight and dephospho- 
rylated by alkaline phosphatase. Internal standards were added and 
the samples were purified by lAC based on a wide range of specific 
monoclonal antibodies against Cks (Faiss et al, 1997). Cks in these 
samples were quantified by UPLC-MS/MS as described above. 

Computational and phylogenetic analysis of the IPT family 

Members of the IPT family in the genomes of 12 selected sequenced 
green plants and algae (Arabidopsis thaliana, Chlamydomonas rein- 
hardtii, Chlorella variabilis, Cyanidioschyzon merolae, Micromonas 
piisilla, Micromonas sp., Oryza sativa, Ostreococcus lucimarinus, 
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Ostreococcus tauri, Selaginella moellendorffii, Volvox carteri, and 
Zea mays) were identified using protein clustering as described 
previously (Bartels et ciL, 2010; Lang et ciL, 2010). The result- 
ing candidate list was collated with the clusters in the Pliytozonie 
plant gene family database (Goodstein et al., 2012) and previous 
descriptions of the IPX family in Arabidopsis, rice, and maize. To 
provide further taxonomic coverage, additional members of the 
IPX family were added based on their InterPro protein domain 
annotation (Hunter et al, 2012). Using the InterPro database web 
interface, proteins containing the IPX protein domain (PF01745; 
IPR002627) from selected cyanobacteria, Alveolates, Euglenozoa, 
Stramenopiles, and Agrobacteriuin tumefaciens were included into 
the candidate set. Sequence data are available as a supplementary 
file at JXB online. Xo derive a high-quality alignment, the manually 
curated alignment of the XRIXl tRNA isopentenyltransferase fam- 
ily from the XreeFam database (XF3 15069) was used as a starting 
point for profile alignment using X-Coffee 8.14 (Notredame et al, 
2000). Based on the resulting full alignment (5517 amino acid col- 
umns), a Neighbor-Joining (NJ) guide tree was calculated using the 
Scoredist matrix (Sonnhammer and HoUich, 2005) implemented 
in a modified version of QuickXree (Howe et al, 2002) applying 
bootstrap resampling with 1000 replicates. Xhe resulting tree was 
used to reorder the alignment and manually curate and clip it using 
Jalview (Clamp et al, 2004) by removing phylogenetically unin- 
formative sites and reducing the full alignment to a clipped and 
reduced ahgnment spanning 203 residues and comprising 97 pro- 
teins. A full list of all IPX family members in this curated alignment 
is provided in Supplementary Xable SI at JXB online. Xhe resulting 
clipped protein alignment was used to infer the IPT gene phytog- 
enies using QuickXree and MrBayes (Ronquist and Huelsenbeck, 
2003). MrBayes was run with parameter settings Iset nst=6 
rates=invgamma Ngammacat=4; prset aamodelpr=mixed; mcmc 
stopval=0.01 stoprule=yes nruns=2 nchains=128 ngen=5000000 
printfreq=1000 samplefreq=100 temp=0.2 starttree^random 
relburnin=yes burninfrac=0.25 savebrlens=yes. 



Results 

Phylogenetic analyses reveal a complex evolutionary 
history of IPTs 

In order to resolve the evolutionary past of the moss IPTs 
and to test their tRNA-IPT ancestry, a comprehensive phylo- 
genetic study was performed considering multiple databases 
and genomes that comprises protein sequences from 22 bac- 
terial (Actinobacteria, Cyanobacteria, and Proteobacteria) 
and 23 eukaryotic species (Amoebozoa, Fungi, Metazoa, 
Rhodophyta, Stramenopiles, and Viridiplantae). The result- 
ing taxonomic distribution of IPT functionality across the 
different kingdoms of Hfe (Fig. 2) raises the question of their 
evolutionary history. 

As previously reported (e.g. Frebort et al, 2011), tRNA- 
IPTs are found in all major clades of bacteria and eukaryotes, 
but are absent from Archaea. In contrast, adenylate-IPTs 
follow a more fragmented distribution. The functionally 
confirmed adenylate-IPTs are found in flowering plants 
(Kakimoto, 2001; Miyawaki et al, 2006; Sakamoto et al., 
2006; Brugiere et al., 2008), the slime mould Dictyostelhim 
discoideum (Taya et al, 1978), and several bacterial plant 
pathogens or symbionts (e.g. Akiyoshi et al, 1985). The pre- 
sent phylogenetic analysis suggests a complex evolutionary 
history of IPTs. Figure 2 depicts the relationships of IPTs 
as a midpoint-rooted phylogram with branch support values 



and colour-coding of branches based on proposed functional 
IPT classes. Besides the tRNA-IPTs, whose function has been 
initially identified as prenylating certain tRNAs in order to 
increase translational precision (Konevega et al, 2006), the 
flowering plant adenylate-IPTs are thought to use ADP or 
ATP as substrate, while the slime mould and bacterial forms 
seem to prefer AMP. Accordingly, Kakimoto (2003) classified 
IPTs into functional groups of DMAPP:tRNA-IPTs (here 
tRNA-IPT, Fig. 2; blue branches), DMAPP:AMP-IPTs (here 
AMP-IPTs, red branches), and DMAPP:ADP/ATP-IPTs 
(here ADT/ATP-IPTs, green branches). 

The absence of adenylate-IPT homologues in Physcomitrella 
gave rise to the hypothesis that in early divergent land plants 
Ck biosynthesis is mediated exclusively by tRNA prenylation 
and subsequent degradation. A functional characterization 
of Physcomitrella IPTl was performed in order to test this 
hypothesis. 

IPT1 is iocaiized in chioropiasts 

Based on sequence analyses, the Physcomitrella IPTl protein 
was predicted to be targeted to chioropiasts (Yevdakova and 
von Schwartzenberg, 2007). For experimental validation of 
the chloroplast localization of IPTl, a IPT1::GFP fusion 
protein was transiently expressed in Physcomitrella proto- 
plasts. The subsequent CLSM-based imaging is consistent 
with a subcellular localization of IPTl in plastids (Fig. 3), 
whereas the control transfection leading to an expression of 
GFP alone showed fluorescence signals in the cytoplasm. It 
was observed that the IPTl ::GFP -mediated fluorescence was 
not evenly distributed within the chloroplast. The confine- 
ment of the GFP signal to the choloroplast has been con- 
firmed by spinning disc CLSM for individual chioropiasts 
(Supplementary Fig. S2 at JXB online). 

Cliaracterization ofd\ipt1 t<noct<out mutants 

In order to draw conclusions on the function of IPTl and 
its relevance for Ck biosynthesis, targeted knockout mutants, 
d|iptl plants, were generated taking advantage of the high 
frequency of homologous recombination in Physcomitrella. 
After selection of stable recombinant plants, the haploid sta- 
tus of the selected strains (#9, #10, and #13) was confirmed 
by flow cytometry (data not shown) in order to exclude that 
protoplast fusion had occurred (Schween et al, 2003, 2005). 
PCR on genomic DNA from the d|iptl plants and subsequent 
sequencing showed that the expected gene replacement, 
namely the partial iptl deletion and full integration of the 
35S::nptII resistance cassette, had occurred at the iptl locus 
(Supplementary Fig. S3a at JXB online). Further, RT-PCR- 
based analyses confirmed the absence of iptl transcript in 
the selected mutants. In order to strengthen the results and 
to minimize unlikely efi"ects of eventual additional random 
ectopic integrations, analyses were performed with three 
independent d|iptl plants; the results for one of the mutants 
(d|iptl#9) is presented in the following. The results for the 
two additional independent mutants confirming the findings 
are given in the supplementary data available at JXB onhne. 
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Fig. 2. IPT phylogeny. Branch lengths and support in the form of posterior probabilities (presented at internal nodes and branch thickness) were Inferred 
using MrBayes (Ronqulst and Huelsenbeck, 2003) based on a manually curated T-Coffee (Notredame et al. , 2000) protein alignment comprising 97 
representative members of the IPT family Confirmed IPT functions are designated with yellow dots: DMAPP:tRNA-IPT references, MiaA_Ecoll (Calllet and 
Droogmans, 1 988), TRIT1 _HUMAN (Golovko et al. , 2000), MOD5„YEAST (Dihanlch ef al. , 1 987), AtlPT2 (Golovko et al. , 2002), AtlPT9 (MIyawakI et al. , 
2006), PplPTI (Yevdakova and von Schwartzenberg, 2007), PplPT3,5 (Patll and Nicanden 201 3); DMAPP:AMP-IPT references, IPT_DCDI (Taya ef al. , 
1978), IPTZ_AGRT7 (AkIyoshI etal., 1985), IPTZ_AGTT5 (Beaty ef a/., 1986): DMAPP:ADT/ATP-IPT references, AtlPTI ,4 (Kaklmoto, 2001), AtlPT3,5,6,7,8 
(MIyawakI ef al. , 2006: Sakaklbara, 2006), OslPTI ,3 (Sakamoto ef al. , 2006), ZmlPT2 (Bruglere ef al. , 2008). Proposed (homology-based) functionality for 
the clades Is given by the colouring of the branches. Gene Identifier and species names are available in Supplementary Table SI at JXB online. 
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Fig. 3. Subcellular localization of IPT1 . Physcomitrella protoplasts 
transfected with (a) empty GFP vector and (b) vector encoding IPT1 ::GFP 
as a C-terminal fusion product. CLSM imaging was performed 5 d after 
transfection. Left, GFP signal; middle, chloroplast autofluorescence; right, 
merged display. Scale bars=10 |xm. 



Phenotype and cytokinin bioassays 
The d|iptl plants expressed a morphological phenotype 
at the level of protonema. Mutant chloronema and cau- 
lonema cells were smaller compared with the wild type 
and showed a higher degree of branching (Fig. 4a, b). 
Compared with the wild type, d|iptl colonies revealed a 
reduced diameter (Fig. 4e, f). d|iptl mutant lines were not 
affected in bud formation; that is, mutant cultures grown 
on solid, as well as in liquid medium still formed buds at 
approximately the same time point and with frequencies 
comparable with the wild type (Fig. 4c, d). Gametophore 
development of the d|iptl mutant strains was also not sig- 
nificantly altered (Fig. 4f). 

In order to check whether the phenotype observed for 
d|iptl mutant protonema could be reversed by externally 
applied Cks, mutants were treated with different concentra- 
tions of iP and BA (20-500 nM). During 2 weeks of micro- 
scopic observation, the phenotype of reduced colony growth 
persisted and budding was induced on the mutant protonema 
to the same extent as on the wild type (Supplementary Fig. 
S4 at JXB onHne). 

Levels of cytokinins in d\ipt1 mutants 

In order to assess the impact of the deletion of iptl on lev- 
els of free and tRNA-bound Cks, they were quantified in 
media and tissue of liquid cultures after 22 d of cultivation 
by UPLC-MS/MS. 

tRNA-bound cytokinins 

Assuming that in Physcomitrella tRNA is a major source 
of Cks, the amount of Ck ribosides in dephosphorylated 
tRNA hydrolysates was determined. cZR had already been 
reported as the dominant Ck type in tRNA of the wild type 
(Yevdakova et al., 2008), correlating well with the domi- 
nance of free cZ-type Cks in tissue. tRNA extracts of d|iptl 
mutants revealed strongly reduced contents of all analysed 
Ck ribosides (Fig. 5a). tZR and DHZR dropped to levels 
below the limit of detection. The more abundant cZR and 
iPR were reduced to <1% compared with the wild type, thus 



wild type d|ipt1 




Fig. 4. Phenotype of d|ipt1 plants (b, d, f, h) in comparison with the wild 
type (a, c, e, g). (a-d) liquid cultures, 19 d old. (a, b) Protonema (scale 
bar=100 |im); (c, d) developing buds (scale bar=20 |xm). (e-h) Development 
on agar plates (scale bar=2mm, 28 d old); (g, h) treated with lOOnM iP. 

revealing the essential importance of IPTl for tRNA A37 
prenylation. 

Free cytokinins in tissue 

In order to uncover the contribution of tRNA-bound Cks to 
the overall Ck production in Physcomitrella, free Cks were 
determined in tissues of the wild type and d|iptl mutants. 
As the methanolic extraction and the applied purification 
and quantification protocol are specific for low molecular 
weight Cks, interference by tRNA-bound Cks is unlikely. 
UPLC-MS/MS analysis revealed that levels of cZ-type Cks 
were reduced from 537 pmol g ' dry weight (DW) in the wild 
type to 132 pmol g ' DW in d|iptl mutants. Furthermore, 
levels of the less abundant DHZ-type Cks (15 pmol g ' DW 
in the wild type) are reduced to only 5 pmol g ' DW in d|iptl 
plants. 

Surprisingly, a different situation was monitored for the 
levels of iP- and tZ-type Cks. In contrast to cZ- and DHZ- 
type Cks, the levels of free iP- and tZ-type Cks were strongly 
increased in the tissues of d|iptl plants (Fig. 5b). For exam- 
ple, iPR5'MP, which at 56 pmol g ' DW was already the 
most abundant iP metabolite in the wild type, was increased 
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Fig. 5. Average levels of isoprenoid cytokinins in 22-day-old liquid cultures 
determined by UPLC-MS/MS. (a) tRNA-bound Cks; two independent 
liquid cultures were harvested and analysed, (b) Free Cks in tissue; three 
independent liquid cultures were harvested and analysed. The results are 
presented as mean values with SDs. ND, not detectable. Data for each 
individual Ck are available in Supplementary Tables S2 and S4 at JXB 
online. (This figure is available in colour at JXB online.) 

~4-fold to levels up to 233 pmol g ' DW in d|iptl mutants 
(Supplementary Table S2 at JXB onhne). 

Extracellular cytokinins 

Previously, the iP-type Cks were found to be highly abun- 
dant in the culture media of Physcomitrella (Reutter et al. , 
1998; Schulz et al., 2001; von Schwartzenberg et al., 2007). 
In the d|iptl plants, the levels of extracellular iP and 
iPRS'MP were determined to be at least 15-fold higher 
than in the wild type. In the mutants, the concentration of 
iPRS'MP in the medium during 22 d of cultivation accu- 
mulated to a concentration of 3697 pM compared with 
242 pM in the wild type (Supplementary Table S3 at JXB 
online). 

In summary, the results of Ck analyses clearly demon- 
strated that the knockout of iptl does not result in a general 



(Lang et al., 2008; Prigge and Bezanilla, 2010), two striking 
aspects are of particular interest in studying Ck biosynthesis 
in the early divergent land plant P. patens: (i) it is the most 
basal member of the green lineage with a sequenced genome 
that encodes major elements of Ck biosynthesis, metabo- 
lism, degradation, and signalling (Rensing et al, 2008; Pils 
and Heyl, 2009); and (ii) the moss genes encoding the first 
step in Ck biosynthesis exclusively represent homologues to 
tRNA-IPTs, a type of IPTs that in flowering plants is consid- 
ered to be only of minor importance to biosynthesis of active 
free Cks (Miyawaki et al, 2006). The origin of these proteins 
and especially the function of IPTl are the key points of 
this study 

Physcomitrella has a unique position In the IPT 
phyiogeny 

Given the topology of the tree and the taxonomic distribu- 
tion generated here (Fig. 2), it can be concluded that the 
origin of the IPT gene family lies within the ancestor of 
extant bacteria. Although the longest internal branch leads 
to AMP-IPTs, this clade is unlikely to represent the ancestral 
form of IPTs, because the AMP-IPT clade is formed entirely 
by proteins from either slime moulds or cyanobacterial and 
proteobacterial plant symbionts or pathogens (red clade in 
Fig. 2). These bacterial genes, encoded by symbiosis or path- 
ogenicity islands (PAIs) within the main genome or plasmids, 
being important for host interactions, are known to have been 
transferred laterally across taxa and consequently have devi- 
ating evolutionary rates (Boyd et al., 2009). This holds true 
also for the putative cyanobacterial adenylate-IPTs as they 
are found only in species that are also observed as symbionts 
of specific fungi, liverworts, hornworts, ferns, gymnosperms, 
and angiosperms (Yamada et al. ,2012, and references within). 
Lastly, the Dictyostelium adenlyate-IPT has been identified as 
a xenologue which was probably acquired by lateral trans- 
fer from bacteria (Eichinger et al, 2005), which is consistent 
with its position in the tree generated here. The true origin 
of bacterial AMP-IPT genes is speculative and could either 
be xenologous or date back to a duplication event in the first 
bacterial Hneage that chose plants as hosts and was subse- 
quently transferred horizontally as part of PAIs. 

Consequently, the origin of IPTs lies within bacterial 
tRNA-IPTs. Based on the tree topology and the taxonomic 
composition of the two plant tRNA-IPT clades (Fig. 2 
class I and II) following the outgroup of AMP-IPTs, the root 
of eukaryotic IPTs most probably traces back to ancestral 
a-proteobacteria and to the initial endosymbiotic event lead- 
ing to extant mitochondria. 

Subsequently, this ancestral mitochondrial protein acquired 
additional domains enabling its function to be extended to 
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the cytoplasm. While Metazoa and Fungi have kept a single 
copy, Plantae, namely the lineage formed by the second endo- 
symbiotic event, either transferred the cyanobacterial copy to 
the nuclear genome or duphcated the mitochondrial tRNA- 
IPT. The topology of the tree favours the latter scenario. The 
ancestral eukaryotic IPT was duplicated and resulted in two 
classes of tRNA-IPT (class I and II). While prasinophyte 
algae [here: Micromonas (Micp) and Ostreococcus (Ostlu 
and Ostta)] as well as seed plants (as exemplified by the well 
characterized Arabidopsis, rice, and maize IPTs) retained 
both classes of tRNA-IPTs, the distinct lineages leading to 
the studied Stramenopiles, the red alga, chlorophyte algae, 
the bryopyhte, and the lycophyte, retained only class I IPTs. 
After an additional duplication event possibly dating back to 
the LCA of flowering plants, one additional copy of class II 
tRNA-IPTs subsequently lost the capability to bind tRNAs 
and diversified into the extant ADP/ATP-IPTs found in flow- 
ering plants (green clade in Fig. 2). 

As previously reported (Sakakibara, 2006; Frebort et al., 

2011) , the seven Physcomitrella IPTs cluster within the 
clade of class I tRNA-IPTs. All other species covered in this 
class code for a single copy tRNA-IPT, which underlines the 
remarkable evolutionary position of Physcomitrella (and 
probably other mosses) regarding IPTs. 

Evidence for a plastidic origin of tRNA-bound Cks 

The locahzation of IPT1::GFP (Fig. 3) corroborated previ- 
ous sequence-based predictions which also suggested plas- 
tidic targeting of IPTl (Yevdakova and von Schwartzenberg, 
2007). Nevertheless, the uneven distribution of the plastidic 
fluorescence signal could also point to a specific suborganel- 
lar localization of the IPTl protein in the stroma, for exam- 
ple in close proximity to the plastid ribosomes (Newell et al. , 

2012) . 

Although the obtained merged pictures and 3D recon- 
structions of CLSM stacks favour the assumption that IPTl 
is present in the stroma of chloroplasts (Supplementary Fig. 
S2 at JXB online), it is not possible at present state to resolve 
its exact subplastidic localization. The observed aggregation 
of IPT1::GFP protein could also be interpreted as an overex- 
pression artefact of the maize ubiquitin promoter resulting 
in an excess of GFP fusion protein which might have precipi- 
tated unevenly in the chloroplast. 

IPTl localized to the chloroplast taken together with 
IPTl as the exclusive mechanism of tRNA prenylation does 
not explain how the tRNA pools in the mitochondria and 
cytosol are prenylated. The present data do not rule out a 
function for IPTl in the cytosol and mitochondria, although 
no GFP signal was detectable in these compartments. An 
inaccurate folding of the fusion protein, for example, could 
have masked targeting sites. For the tRNA-IPT MODS from 
Saccharomyces cerevisiae, it has been shown that according to 
different translational starts, the protein is either transported 
into mitochondria or stays in the cytosol (Gillman et al., 
1991). Further, it has been shown for AtIPT3 in Arabidopsis 
that farnesylation can also determine targeting to different 
subcellular compartments (Gahchet et al, 2008). Due to the 



scarce knowledge of tRNA trafficking in plant cells, it can 
only be speculated on as to whether A37 modified tRNAs 
can be translocated from the chloroplast to other cellular 
compartments. 

Under the assumption that the targeting of IPT1::GFP 
was not biased by artefacts, the present results point to the 
chloroplast providing a large amount of Physcomitrella pre- 
nylated tRNA. This is a contrast to AtIPT2, a tRNA-IPT 
from Arabidopsis, which has been shown to be localized in 
the cytosol where it is considered to be a significant source 
of cZ-type Cks using DMAPP from the mevalonate pathway 
as a side chain donor (Kasahara et al, 2004). For the other 
tRNA-IPT of Arabidopsis, AtIPT9, to which PpIPTl is an 
orthologue, no information on localization is available. 

iPTI is tlie main enzyme catalysing tRNA prenyiation 

The deletion of iptl as one of seven ipt genes had a severe 
effect on the prenylation of A37 in tRNA. The amount of 
tRNA-bound cZR, iPR, tZR, and DHZR in the d|iptl 
mutants was strongly reduced (Fig. 5) and thus the contribu- 
tion of the remaining six IPTs to tRNA prenylation can be 
considered as insignificant in moss protonema. It is therefore 
concluded that plastidic IPTl is almost solely responsible for 
this A37 modification in Physcomitrella tRNA. The obtained 
results are comparable with the report of Dihanich et al. 
(1987) who found in yeast that the deletion of the single copy 
tRNA-IPT gene Mod5 led to residual levels of <1.5% tRNA- 
bound iPR compared with the wild type. 

If the tRNA pathway was the only pathway for Ck bio- 
synthesis, as suggested by the apparent absence of adenylate- 
IPTs, a Ck-deficient tRNA as obtained after IPTl knockout 
would have led to an overall Ck-deficient plant. These plants 
should express strongly impaired bud formation as shown for 
Ck deficiency after overexpression of cytokinin oxidase (von 
Schwartzenberg et al. , 2007). 

Indeed the IPTl knockout plants showed a strong phe- 
notype which was mainly characterized by altered colony 
growth, cell shape, and branching of protonema (Fig. 4a, b). 
This phenotype could not be reversed by addition of exog- 
enous Cks (Fig. 4; Supplementary Fig. S4 at JXB online). 
However, the bud formation as a sensitive marker of Ck 
response in Physcomitrella was neither delayed nor reduced in 
comparison with the wild type. The reduced growth in colony 
area of the mutants with its concomitantly enhanced differ- 
entiation (Fig. 4e, f) can be mimicked in wild-type agar cul- 
tures by exogenous application of iP (Fig. 4g, h) (Thelander 
et al., 2005; Richter et al., 2012). This morphological feature 
of d|iptl mutants indicates an overproduction rather than the 
expected limitation in biosynthesis of bud-inducing Cks. 

Profiling of free C/cs indicates a tRNA- independent Cl< 
biosynthesis 

A fully satisfactory explanation of the d|iptl phenotype is 
only possible after a detailed analysis of the profile of free 
Cks in the mutant and wild type, which uncovered a differ- 
ential situation with respect to individual Ck types. While 
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the levels of tRNA-bound Cks (cZR, iPR, tZR, and DHZR) 
were strongly reduced, the level of the free Cks revealed a sig- 
nificant reduction only for cZ- and DHZ-type Cks to about a 
quarter of the respective amount in the wild type. Since IPTl 
is localized in chloroplasts, it can be deduced that the major- 
ity of free cZ- and DHZ-type Cks is of plastidic tRNA origin. 

In contrast to the reduction of free cZ- and DHZ-type Cks, 
surprisingly a 4-fold increase for the iP- and tZ-type Cks was 
found. The fact that there is no overall loss of free Cks in the 
mutants but an increase of tZ and iP conjugates coupled to 
an almost total loss of tRNA-bound Cks, is strong evidence 
for a an unexpected tRNA-independent pathway for Ck bio- 
synthesis in Physcomitrella. 

This finding is in obvious disagreement with the predicted 
ipt gene complement, which seems to be comprised exclusively 
of homologues of tRNA-IPTs. Thus, the biosynthetic origin 
of tZ- and iP-type Cks needs to be re-assessed as tRNA is 
very unlikely to be their only source. 

In Arabidopsis it was described that different pathways 
are responsible for the formation of distinct types of Cks. 
The tRNA-IPT-deficient double mutants of Atipt2,9 were 
shown, like the d|iptl mutant, to contain no Cks in their 
tRNA fraction (Miyawaki et al, 2006). However, in contrast 
to the Physcomitrella mutants, with 25% residual free cZ- 
type Cks, the AtIPT2,9-deficient plants did not produce any 
cZ-type Cks. 

In contrast to Arabidopsis where tZ- and iP-type Cks are 
dominant, the cZ type are the most abundant cytokinins in 
Physcomitrella. Despite this difference in the Ck profiles, 
in both plants iP and tZ represent the physiologically most 
active forms (Sakakibara, 2006; von Schwartzenberg et al, 
2007). 

Physcomitrella and Arabidopsis differ greatly in the levels 
of free iP- and tZ-type Cks in the absence of one (PpIPTl) 
or two tRNA-IPTs (AtIPT2 and -9). While in Arabidopsis 
the level of free iP- and tZ-type Cks is not greatly affected 
in the tRNA-IPT double mutants, the d|iptl mutants show 
a vastly increased level of those Cks, thus hinting at dif- 
ferences in Ck biosynthesis, and its regulation, between 
the bryophyte Physcomitrella and flowering plants. In 
Physcomitrella, the deficiency of IPTl apparently causes 
an up-regulation of the tRNA-independent formation of 
iP and tZ. 

The knockout of iptl seems to interfere mainly with the bio- 
synthesis but not with the interconversion between different 
conjugates of free Cks, thus leading to either generally increased 
or reduced overall concentrations of the different Ck types 
(Fig. 5b; Supplementary Table S2 at JXB online). The relative 
distribution of the different Ck forms within one type (e.g. the 
ratio of bases to ribosides or nucleotides) is only slightly affected 
(Supplementary Fig. S5, Table S2 available at JXB online). 

Concerning the profile of free Cks, the deletion of IPTl 
leads to a shift from cZ-type dominance in the wild type 
to an iP-type dominance of free cytokinins in d|iptl tissue 
(Fig. 5b). These strong changes in the internal Ck profile also 
affect the concentration of Cks in the culture medium, with 
the concentration of iP-type Cks being greatly enhanced in 
the mutants. 



The increase of iP-type Cks in the culture medium is likely to 
explain certain aspects of the phenotype of d|iptl plants. Both 
a reduction in colony growth area and an increased differen- 
tiation might be a direct consequence of the increased external 
levels of iP and iPR, mirrored by the growth habitus of wild- 
type plants grown on iP-containing media. The fact that the 
frequency of bud induction is comparable in the mutant and 
wild type (Fig. 2; Supplementary Fig. S4b at JXB online) shows 
that the d|iptl plants are probably not affected in Ck signalling. 

The results of the analyses of d|iptl mutants suggest that in 
Physcomitrella tRNA-bound cZR is most probably the main 
source for free cZ-type Cks. Considering the loss of tRNA- 
bound Cks and the concomitant existence of substantial 
amounts of free iP- and tZ-type Cks, the IPTl -deficient plants 
strongly indicate the existence of a second, tRNA-independ- 
ent pathway in the moss. This pathway seems to be important 
in the production of the active Cks in Physcomitrella. 

Given the Ck profile of the d|iptl mutant and the hence 
deduced existence of a tRNA-independent Ck biosynthetic 
pathway, the clustering of moss IPTs separate from adenylate- 
IPTs leads to three possible hypotheses, (i) There is an alter- 
native pathway mediating the tRNA-independent production 
of Cks by an as yet unknown enzyme, (ii) One or several of 
the other six moss IPTs belong to an as yet uncovered family 
of (maybe bryophyte-specific) IPTs which have convergently 
lost the capacity to bind tRNA, evolving adenylate-IPT func- 
tionality in parallel to flowering plant ADP/ATP-IPTs. (iii) 
One or several of the moss IPTs indeed belongs to the same 
ADP/ATP-IPT family as those of flowering plants, but the 
phylogenetic signal connecting them has been obscured by an 
as yet undetermined mechanism. In the last case, the origin 
of plant adenlyate-IPTs would trace back to the ancestor of 
land plants. Based on the available data, especially the topol- 
ogy and branch lengths of the PpIPT subtree, alternative (ii) 
seems to be the most probable evolutionary scenario. 

The functionally confirmed tRNA-IPT PpIPTl clusters 
basal to all other moss IPTs. Some of these six additional 
IPTs deviate substantially in their branch lengths, indicat- 
ing changes in the rate of evolution. In order to see whether 
the increased amounts of iP- and tZ-type Cks are due to an 
increased expression of one of the remaining IPTs in the 
d|iptl mutant, real-time analyses have been performed, iptl.l 
and ipt2.2 did not show expression either in the wild type or 
in the mutant. The remaining ipt genes {ipt3-ipt6) showed 
no strong or consistent changes (Supplementary Fig. S6 at 
JXB online). Which of the proposed evolutionary hypotheses 
ultimately can be verified clearly requires additional experi- 
mental work. Therefore, studies have been initiated to char- 
acterize functionally the in planta function of the remaining 
members of the Physcomitrella IPT family in order to clarify 
the identity of the unexpected tRNA-independent Ck bio- 
synthetic pathway observed in Physcomitrella. 



Supplementary data 

Supplementary data are available at JXB online. 
Figure SI. Vector card for pLNU-GFP_Ppiptl. 
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Figure S2. Localization of IPTl within Physcomitrella 
chloroplasts. 

Figure S3. Generation and characterization of d|iptl 
mutants. 

Figure S4. Ck response of d|iptl mutants and the wild type. 
Figure S5. Relative level of cytokinins in 22-day-old Hquid 
cultures. 

Figure S6. Relative expression of the ipt gene family in 
d|iptl mutants. 

Table SI. Gene identifiers for the phylogenetic analyses 
performed. 

Table S2. Average levels of intracellular free isoprene-type 
Cks in tissue. 

Table S3. Average levels of extracellular free isoprene-type 
Cks in media. 

Table S4. Average levels of tRNA-bound isoprene-Ck ribo- 
sides in tissue. 
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